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Abstract 

An updated version of the PNJL model is used to study the thermodynam- 
ics of Nf = 2 quark flavours interacting through chiral four-point couplings 
and propagating in a homogeneous Polyakov loop background. Previous PNJL 
calculations are extended by introducing explicit diquark degrees of freedom 
and an improved effective potential for the Polyakov loop field. The mean 
field equations are treated under the aspect of accommodating group theoret- 
ical constraints and issues arising from the fermion sign problem. The input 
is fixed exclusively by selected pure-gauge lattice QCD results and by pion 
properties in vacuum. The resulting (T, fx) phase diagram is studied with 
special emphasis on the critical point, its dependence on the quark mass and 
on Polyakov loop dynamics. We present successful comparisons with lattice 
QCD thermodynamics expanded to finite chemical potential fi. 



1 Introduction 



Reconstructing the phase diagram and thermodynamics of QCD in terms of field 
theoretical quasiparticle models is an effort worth pursuing in order to interpret 
lattice QCD results [1-7] and extrapolate into regions not accessible by lattice com- 
putations. A promising ansatz of this sort is the PNJL model [8—11], a synthesis 
of Polyakov loop dynamics with the Nambu & Jona-Lasinio model, combining the 
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two principal non-perturbative features of low-energy QCD: confinement and spon- 
taneous chiral symmetry breaking. This paper extends our previous PNJL calcula- 
tions [10, 11] in several directions. First, diquark degrees of freedom are explicitly 
included. Diquark condensation at large quark chemical potential is explored in the 
presence of a Polyakov loop background. Secondly, in comparison with our previous 
work, the effective potential which controls the thermodynamics of the Polyakov 
loop field is improved such that group theoretical constraints are implemented, fol- 
lowing Ref . [8] . Issues of the mean field approximation in the context of the fermion 
sign problem are discussed in comparison with previous work. 

The aim of the present paper is to investigate the phase diagram resulting from 
this approach. Of special interest is the location of the critical point, its dependence 
on the quark mass and the role of the Polyakov loop as indicator of the deconfinement 
transition. Predictions for the leading coefficients in a Taylor expansion of the 
pressure in powers of the quark chemical potential will turn out to be remarkably 
successful in comparison with corresponding lattice QCD results. 



2 The PNJL model 

The two-flavour PNJL model (now including diquark degrees of freedom) is specified 
by the Euclidean action 

r p=\/T r 

S E (*P, ^, 0) = J drj d 3 x 8 T i, + H{$, 4>\ 0)] + SSe^, T) (1) 

with the fermionic Hamiltonian density 1 : 

H = -#t(a- V + 7 4 m o -0)^ + V(^,^ t ) , (2) 

where ip is the Nf — 2 doublet quark field and m = diag(m u , ma) is the quark mass 
matrix. The quarks move in a background color gauge field = At = iA , where 
A = 5^ gA£t a with the SU(3) C gauge fields A£ and the generators t a = \ a /2. The 
matrix valued, constant field relates to the (traced) Polyakov loop as follows: 



Pexp(i drA^j = ^Tre* /T , (3) 



In a convenient gauge (the so-called Polyakov gauge), the matrix is given a diagonal 
representation 

= 03 A 3 + 08 A 8 , (4) 

which leaves only two independent variables, 03 and 8 . The piece 5Se = \U 
of the action (1) controls the thermodynamics of the Polyakov loop. It will be 



1 a = 7o 7 and 74 = 270 in terms of the standard Dirac 7 matrices. 
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specified later in terms of the effective potential, W($,T), determined such that 
the thermodynamics of pure gauge lattice QCD is reproduced for T up to about 
twice the critical temperature for deconfinement. At much higher temperatures 
where transverse gluons begin to dominate, the PNJL model is not supposed to be 
applicable. 

The interaction V in Eq. (2) includes chiral SU(2) x SU(2) invariant four- 
point couplings of the quarks acting in pseudoscalar-isovector/scalar-isoscalar quark- 
antiquark and scalar diquark channels: 



-4 



H 



(W) +$i-y 5 T'i/>) --[{^C l5 T 2 X 2 ^ T )(Vl 5 r 2 X 2 C^)} , (5) 



where C is the charge conjugation operator. We can think of Eq. (5) as a subset 
in the chain of terms generated by Fierz-transforming a local color current-current 
interaction between quarks, 

In this case the coupling strengths in the quark-antiquark and diquark sectors are 
related by G = the choice we adopt. The minimal ansatz (5) for V is motivated 
by the fact that spontaneous chiral symmetry breaking is driven by the first term 
while the second term induces diquark condensation at sufficiently large chemical 
potential of the quarks. Additional pieces representing vector and axialvector qq 
excitations as well as color-octet diquark and qq modes are omitted here. We have 
checked that their effects are not important in the present context. 

The NJL part of the model involves three parameters: the quark mass which 
we take equal for u- and c?-quarks, the coupling strength G and a three-momentum 
cutoff A. We take those from Ref.[10]: 

4 

m u4 = 5.5 MeV , G = -H = 10.1 GeV" 2 , A = 0.65 GeV , 

which were fixed to reproduce the pion mass and decay constant in vacuum and the 
chiral condensate as m w = 139.3 MeV, f„ = 92.3 MeV and (^„^ u ) = -(251 MeV) 3 . 

The effective potential W(<3>, T) which controls the dynamics of the Polyakov loop 
has the following properties. It must satisfy the Z(3) center symmetry of the pure 
gauge QCD Lagrangian. In the low-temperature (confinement) phase W($,T) has 
an absolute minimum at $ = 0. Above the critical temperature for deconfinement 
(T ~ 270 MeV according to pure gauge lattice QCD results) the Z(3) symmetry is 
spontaneously broken and the minimum of W($,T) is shifted to a finite value of $. 
In the limit T — > oo we have $ — > I. 

In our previous Ref. [10] the simplest possible polynomial form was chosen for U. 
In the present work an improved expression, guided by Ref. [8], replaces the higher 
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order polynomial terms in $, $* by the logarithm of the Jacobi determinant 

which results from integrating out six non-diagonal Lie algebra directions while 
keeping the two diagonal ones, 03,8, to represent <3>. This suggests the following 
ansatz for U : 

U ^ T P = --a(T) <F$ + b(T) In [l - 6 $*$ + 4 (<T 3 + $ 3 ) - 3 ($*$) 2 ] (6) 
with 

a(T) =a + ai (^j + a 2 (^j ' , b(T) = b 3 (^j ' . (7) 

With its logarithmic divergence as $, $* — > 1, this ansatz automatically limits the 
Polyakov loop $ to be always smaller than 1, reaching this value asymptotically only 
as T — > oo. Following the procedure as in [10], a precision fit of the parameters Oj 
and 63 is performed in order to reproduce lattice data for pure gauge QCD thermo- 
dynamics and for the behaviour of the Polyakov loop as a function of temperature. 
The critical temperature T for deconfinement in the pure gauge sector is fixed at 
270 MeV in agreement with lattice results. 

The results of this combined fit are shown in Fig. 1 and the dotted line of Fig. 2. 
The corresponding parameters are 

a = 3.51 , ai = -2.47 , a 2 = 15.2 , 63 = -1.75 . 

The fit was constrained by demanding that the Stefan-Boltzmann limit is reached 
within the model at T — > 00 and by enforcing a first-order phase transition at 
To. The first constraint determines ao = ^§|- ~ 3.51. The second constraint fixes 
b 3 = —0.108 (a + ai + a 2 ). The two remaining parameters a± and a 2 were determined 
using a mean square fit. In this fit the Polyakov loop data set [6] was given a stronger 
weight than the pressure, energy density and entropy [5]. This was done in order 
to counterbalance the smaller number of Polyakov loop data points against the 
larger number of data sets for pressure, energy density and entropy. The resulting 
uncertainties are estimated to be about 6% for a x and less than 3% for a 2 . These 
independent errors propagate to 63 giving an uncertainty of about 2 %. 

Next, the PNJL action is bosonized and rewritten in terms of the auxiliary 
scalar and pseudoscalar fields (a, 7?), and diquark and antidiquark fields (A, A*). 
The thermodynamic potential of the model is evaluated as follows: 

where the sum is taken over Matsubara frequencies u n = (2n + l)irT. The inverse 
Nambu-Gor'kov propagator including diquarks is: 

5-i a, , tt\ - ( *7o u n ~ 7 • P - m + 7o - i<i>) A7 5 r 2 A 2 

^ -A 7 5 t 2 A 2 0J n - 7 • p - m - 70 {/i - i<p) 

(9) 
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Figure 1: Left: Fit to scaled pressure, entropy density and energy density as functions of the 
temperature in the pure gauge sector, compared to the corresponding lattice data taken from 
Ref. [5]. Right: Resulting effective potential (6) that drives spontaneous Z(3) symmetry breakdown 
at T = To- 



Just as in the standard NJL model, quarks develop a dynamical (constituent) mass 
through their interaction with the chiral condensate: 

m = m — (a) = m — G{ipip). (10) 

With the input parameters previously specified one finds m = 325 MeV at T = 0. 

Note that introducing diquarks (and anti-diquarks) as explicit degrees of freedom 
implies off-diagonal pieces in the inverse propagator (9). As a consequence, the 
traced Polyakov loop field $ and its conjugate $* can no longer be factored out 
when performing the Trln = lndet operation in the thermodynamic potential (8), 
unlike the simpler case treated in our previous Ref. [10]. The explicit evaluation of 
energy eigenvalues now involves 03 and 8 as independent fields. The final result 
for Q is then given as: 

a 2 A* A 

a = u(*,T) + - + — 

~ 2N ' I (Sp£{ rln[l + e ~ £ ' /T] + 1 2 AE '} ■ (U) 

The difference 0,(T) — Q(T = 0) is to be used in the actual calculations. The quasi- 
particle energies Ej, denoted by indices j running from 1 to 6, have the following 
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explicit expressions with e(p) = a/p 2 + m2 '- 



Ei, 2 = e(p)Th, 

£5,6 = N /(e(p)-/i r ) 2 + |A| 2 =Fi03, (12) 



with 



fl b = fi + 2l -j= , H r = fl-l-f=. (13) 

In Eq.(ll), A£j = £j — £ ± ji is the difference of the quasiparticle energy Ej and 
the energy of free fermions (quarks), where the upper sign applies for fermions and 
the lower sign for antifermions. It is understood that for three- moment a \p\ above 
the cutoff A where NJL interactions are "turned off', the quantities a and A, A* 
are set to zero. 

The thermodynamical potential (11) involves the bosonic field variables a, A, 
03 and 08- In the mean field theory the integration over all field configurations 
{a, A, 03, 08} in the calculation of the partition function is approximated by a single 
field configuration, {a, A, 3 , 0s} m f • For a purely real action the optimal mean field 
configuration is the one which determines a minimum of Q = Q(a, A, 3 , 8 ). The 
necessary condition for this is 

9n =0. 



9 (<T, A, 03, 8 ) 

Here the action is treated in analogy with a Landau effective action which identifies 
the fields with approximate order parameters. In case of the PNJL model, however, 
the thermodynamical potential Q is complex in the presence of the Polyakov loop 
background and at non-zero chemical potential A minimization of a complex 
function as such is void of meaning. This is the so-called fermion sign problem in 
the present context. However, even for a complex Euclidean action Se, one can still 
search for the configuration with the largest weight in the path integral and refer to 
this as the mean field configuration. An analysis of the phase and the absolute value 
of the weight e~ SE immediately shows that this mean field configuration maximizes 
|e~ 5E | and consequently minimizes Ref2. The mean field equations derived from fl 
are 

dReQ 

9(a,A,0 3 ,0 8 ) " ' ( } 

the condition we adopt. 

In previous publications [10,12,13], the last two mean field equations, = 
■§§- = Q, were replaced by |§ = J|p- = 0. These relations are in principle equivalent, 
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given that $ = $(03, 0s) and $* = $*(03, 0s) a s implied by Eq. (3). The constraint 
under which such a change of variables can be done is that the temporal gauge fields 
remain real quantities: 03, 0s G R. 

Abandoning this constraint would introduce different chemical potentials for 
quarks of different colors. 2 Using equation (3) it can easily be derived that in 
the case where 03 G K and 8 G M, $ and <£>* genuinely have to be the complex 
conjugate of each other. 

The (thermal) expectation values (<&) and (<&*) of the conjugate Polyakov loop 
fields must be real quantities as argued in [14]. This applies to the PNJL model as 
well, in the sense that the mean field action iSmf changes into its complex conjugate 
under charge conjugation. Enforcing ($) G K. and (<£>*) G K means $ = $* for the 
mean field configurations that satisfy Eq. (14). With the constraint of 03 and 08 
being real, this implies 08 = leaving only 03 as an independent variable. 

In previous work [10, 12, 13] ($) and ($*) have been treated as independent 
real quantities in the minimization of Q. This procedure, without the constraints 
imposed by 03, 08 G R, tends to overestimate the difference between ($) and ($*). 
While this unphysical feature has only marginal consequences for global properties 
such as the pressure, it does have a visible influence on more detailed quantities as 
discussed in section 3.1 

Within the present mean field context defined by Eq. (14), fluctuations beyond 
mean field are at the origin of ($) ^ ($*) for /i ^ 0. This paper deals with self- 
consistent solutions and predictions of the mean- field equations (14). While further 
extensions including quantum fluctuations [15] will be subject of a forthcoming 
publication [16], we can already anticipate one of the results, namely that the effects 
of fluctuations, leading to ($) 7^ ($*) at finite chemical potential, turn out not to be 
of major qualitative importance in determining the phase diagram. This forthcoming 
publication will present a way how to deal with the fermion sign problem in the 
context of the PNJL model. Discussions of the fermion sign problem in other models 
can be found in [14, 17]. 

3 Results 

Solution of the mean-field equations (14) yields the chiral condensate, (ipip) = cr/G, 
the color- antitriplet diquark condensate, A, and the Polyakov loop exponent 03 as 
functions of T, \i. The resulting prediction for the traced Polyakov loop $ at /i — 
(where the diquark condendate vanishes, A = 0) is shown in Fig. 2 (continuous 
line) in comparison with the corresponding lattice data taken from Ref. [7] (full 
symbols). The agreement is quite remarkable. The improvement in comparison to 

2 This can formally be seen when using the simple prescription [i — > [i — i (<f> 3 A3 T + <f> s A 8 T) 
to do the transition from the NJL- to the PNJL-Nambu-Gor'kov propagator. 
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Figure 2: Using the fit of the Polyakov loop (dotted line) to lattice results taken from Ref. [6] 
in the pure gauge sector (empty symbols) , the PN JL model predicts the Polyakov loop behaviour 
as a function of temperature in the presence of dynamical quarks (solid line). This prediction is 
compared to lattice data in two flavours (full symbols) taken from Rcf. [7]. 



previous calculations is primarily due to the improved ansatz for the Polyakov loop 
potential. In the presence of quarks, the deconfinement transition is no longer first- 
order as in pure gauge QCD. It becomes a smooth crossover when quarks couple to 
the Polyakov loop field. The critical temperature for deconfinement is now decreased 
from 270 MeV to a smaller value 3 around 215 MeV (not evident from Fig. 2 where the 
results are plotted as functions of T/T c ). In Fig. 3 we show in addition the predicted 
temperature dependence of the two-flavour chiral condensate (ipip) in comparison 
with lattice data [18]. 

3.1 Finite chemical potential 

Lattice results at finite quark chemical potential are obtained as Taylor expansions of 
the thermodynamical quantities in the parameter \ijT around zero chemical poten- 
tial. Here we perform the same kind of expansion in the PNJL model and compare 
with Taylor coefficients deduced from lattice data. Examples are the coefficients in 
the expansion of the pressure p = —Q: 



3 Note that the critical temperature in full lattice QCD reported in [7] is T c = 202 MeV. 




n=0 



with 




(15) 
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1. 




Figure 3: The spontaneous chiral symmetry breaking mechanism of the PNJL model generates 
a temperature dependent chiral condensate (4>ip) (solid line), which is compared here to lattice 
QCD results in two flavours shown in Ref. [18]. 



and even n. Specifically: 

1 d 2 (p/T 4 ) 



c 2 



c 6 



2 d^/Tf 
1 d 6 (p/T 



c 4 



1 d 4 (p/T 4 



/i=0 



720 d(fi/ry 



c 8 



24 <9(/VT) 4 
1 d 8 (p/T 



fl=0 



40320 d(/i/TY 



(16) 



At=0 



Results for these coefficients are shown in Fig. 4. We notice in particular the remark- 
ably good agreement between the calculated " susceptibility" C4 and the lattice data. 
This quantity has recently been computed in Ref. [12] using the previous version of 
our PNJL model, Ref. [10], with a less satisfactory outcome. Now, with the improved 
effective potential U as described in Eq.(6), the agreement is significantly better. We 
suspect that previous calculations did not approach the Stefan-Boltzmann limit in 
an acceptable way [19], due to the large gap between $ and that persisted up 
to rather high temperatures. As discussed above this large split occurs when not 
properly keeping all constraints on 3 and (fis under control. 



3.2 Phase diagram 

We now turn to the phase diagram in the (T, /i) plane as derived from this updated 
version of the PNJL model. The left panel of Fig. 5 shows the phase diagrams 
in the (T, //)-plane computed using the PNJL model in comparison with the NJL 
model (the limiting case in which $ = 1). Of particular interest is the location of 
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Figure 4: Second, fourth, sixth and eighth moments of the pressure difference with respect to the 
chemical potential, plotted as functions of the temperature. (Note that the temperature scales of 
the upper and lower graphs are different.) We compare to lattice data (diamonds with errorbars) 
taken from Ref. [4] 
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the critical endpoint at which the chiral and deconfinement crossover transitions at 
lower fi turn into a first-order phase transition above some critical \x. The crossover 
is not a phase transition. Therefore there exist several ways to locate the position 
of a crossover transition. In the present calculations the crossover line is determined 
using the order parameters in the chiral limit (the chiral condensate) and the pure 
gauge theory (the Polyakov loop) respectively. Since these order parameters show 
their strongest changes as functions of temperature along the crossover transition 
lines, we determine their position by local maxima of da/dT and d$/dT 4 . 

The crossover transition lines fixed by either the susceptibilities of a and <3> or by 
maximal changes with temperature, i.e. zeros of d 2 a/dT 2 or d 2 $/dT 2 , do coincide 
with the critical point for our PNJL model in the absence of diquarks (see lower 
panel of Fig. 5). This is a consequence of the divergences in these quantities at the 
critical point. However, when including diquarks, a coincidence of critical point and 
crossover transition line is not guaranteed. 

One finds that the critical endpoint depends sensitively on the degrees of freedom 
involved. From its position in the restricted NJL case (see also [20]) this point 
is shifted to higher T by both, the effective Polyakov loop potential, and by the 
presence of diquark degrees of freedom. Near the critical endpoint not including 
diquarks, ^ diverges together with the chiral susceptibility. This extreme behaviour 
is not observed in the case with inclusion of diquarks. The region where this critical 
behaviour would appear is now already located in the diquark dominated phase. 

Thus there is a qualitative difference of the critical endpoints in these two com- 
pared cases: not including diquarks the critical endpoint lies on top of the merging 
chiral and deconfinement crossover transition lines, while in the case including di- 
quarks the critical endpoint is shifted away from this line. The critical endpoint now 
lies on the second order transition line bordering the diquark dominated phase (see 
lower panel of Fig. 5). I. e. the endpoint is not at the junction of all three transition 
lines and therefore is not a tri-critical point but still a critical point. 

Next we use the PNJL model including diquark degrees of freedom to study the 
dependence of the position of the critical endpoint on the bare (current) quark mass. 
The upper right panel of Fig. 5 shows phase diagrams in the chiral limit, for current 
quark masses itlq = 5.5 MeV and itlq = 50 MeV. The change of the critical endpoint 
with varying quark mass mainly reflects the dependence of the critical chemical 
potential on the quark mass. The presence of the diquark dominated phase appears 
to stabilize the temperature of the critical endpoint at rather high values. 

Generally, the PNJL model generates the critical endpoint at a temperature 
which is significantly higher than the one found with the standard NJL model, i. e. 
ignoring Polyakov loop dynamics. The reason is that the diquark phase as well as 

4 Other frequently used and closely related criteria for the definition of crossover transition lines 
involve chiral or Polyakov loop susceptibilities. This does not lead to any significant differences for 
the phase diagram in comparison with the method applied here. 
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Figure 5: Upper left panel: comparison of the phase diagrams of NJL and PNJL model. The 
cross-over of the chiral condensate is drawn solid, first-order lines are dashed and second order lines 
dotted. Upper right panel: comparison of the phase diagram at different current quark masses 
with inclusion of diquark degrees of freedom. (Note the scale on the temperature axis.) Lower 
panel: comparison of the PNJL model with and without inclusion of diquarks. 
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Figure 6: Dependence of the gap A on the presence of the Polyakov loop. The solid lines are the 
solutions to the self consistency equations of the NJL and the PNJL model at /j, = 0.4 GeV. The 
dashed lines are obtained by enforcing fixed values for the Polyakov loop. Note that the PNJL 
model with the Polyakov loop fixed at <f> = 1 (deconfinement) coincides with the self consistent 
solution of the NJL model. 



the chiral phase is stabilized by the confinement emulation via the effective Polyakov 
loop potential. The size of the gap A is strongly influenced by the Polyakov loop. 
The detailed dependence of the gap on the Polyakov loop is displayed in Fig. 6. 
The systematics of this effect becomes evident when the Polyakov loop is held at 
fixed values and varied. The gap resulting from this calculation is then compared 
to the gap in the PNJL model (with self-consistent determination of $) and in the 
NJL model. The case where the Polyakov loop is fixed to $ = 1 (i. e. complete 
deconfinement) coincides with the NJL calculation. 

The presence of the Polyakov loop restricts the phase space available for quarks 
in the vicinity of their Fermi surface where Cooper pair condensation takes place. 
Hence a higher temperature is effectively required to break the pairs. This is the 
primary reason for the difference in behavior of the gap A when comparing NJL and 
PNJL results in Fig. 6. 



3.3 Speed of sound 

The velocity of sound v s , determined by 



1 dtt 



v=- 



C v dT 



;i7) 



v 



with the specific heat Cy = —T(d 2 Q/dT 2 )v, shows a pronounced dip near the chi- 
ral and the deconfinement transition. This local minimum of the speed of sound 
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Figure 7: Left panel: Squared speed of sound (dashed) and ratio | (solid), and their dependence 
on the temperature at vanishing chemical potential. The two panels to the right: Studies of the 
sound velocity in the vicinity of the critical endpoint in the PNJL model without diquarks (at 
chemical potential /i = /icEp)- The PNJL model without diquarks is drawn with solid lines, the 
PNJL model with diquarks with dashed lines. Center panel: square of the velocity of sound. Right 
panel: pressure over energy density |. 



becomes deeper in the vicinity of the critical endpoint of the first-order phase tran- 
sition line, separating the chiral phase at low chemical potential (// < 1.5 T c ) from 
the diquark phase at high chemical potential (fx > 1.5 T c ). When neglecting diquark 
degrees of freedom the speed of sound vanishes at the critical endpoint (solid curve 
in the central panel of Fig. 7). 

Correspondingly, the specific heat diverges at this point. When diquark degrees 
of freedom are included in the calculation the critical endpoint is shifted such that 
the region of vanishing speed of sound would already be placed within the diquark 
dominated phase. This is why a vanishing speed of sound is not observed in the 
model including explicitly diquark degrees freedom (dashed curve in the central panel 
of Fig. 7). The discontinuity at higher temperatures is generated by the second order 
phase transition separating the diquark dominated phase from the high temperature 
phase. Above this transition the two versions of the PNJL model (with and without 
explicit diquarks) become equivalent. 

4 Concluding remarks and outlook 

We have pointed out that an updated version of the PNJL model over and beyond 
the one used in [10, 12] leads to significantly better agreement with lattice data [4, 7], 
especially when extrapolating to finite chemical potential \i. The combination of only 
two principal ingredients: chiral symmetry restoration and an effective potential 



14 



ansatz for the confinement order parameter, appears to be sufficient to reproduce 
the available full QCD lattice computations to an astonishingly high accuracy, at 
least for temperatures T up to about 2T C . The improvements shown in this paper 
in comparison with previous results [10, 12] originate in a better representation of 
the Polyakov loop part of the PNJL model. Taking into account the proper SU(3) 
constraints is crucial for an effective description of the thermodynamical implications 
of confinement. 

Incorporating explicit diquark degrees of freedom influences the position and the 
nature of the critical endpoint in the (T, fj,) phase diagram. The critical endpoint in 
the presence of diquarks is the connecting point between the chiral crossover tran- 
sition line and the second order transition bordering the diquark dominated phase, 
while in the absence of diquarks it is the junction point of the chiral and decon- 
finement crossover transition. The critical point in the PNJL model with diquarks 
turns out not to coincide with the critical (diverging) behaviour of susceptibilities 
related to the chiral condensate and the Polyakov loop. 

Further developments now aim for an extension of the present framework to Nf = 
3 in order to explore the rich structure of colour superconducting (diquark) phases 
with three quark flavours and the additional effects of Polyakov loop dynamics. 
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